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Abstract. SGR 1806-20 has been observed for more than 2 years with the INTEGRAL satellite. In this period 
the source went from a quiescent state into a very active one culminating in a giant flare on December 27 2004. 
Here we report on the properties of all the short bursts detected with INTEGRAL before the giant flare. We 
derive their number-intensity distribution and confirm the hardness-intensity correlation for the bursts found by 
|Gotz et al. (2004a) | Our sample includes a very bright outburst that occurred on October 5 2004, during which 
over one hundred bursts were emitted in 10 minutes, involving an energy release of 3xl0 42 erg. We present a 
detailed analysis of it and discuss our results in the framework of the magnetar model. 
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1. Introduction 

SGR 1806-20 is currently one of the most active Soft 



Gamma-ray Repeaters. These sources (see Hurley (2000) 



Woods & Thompson (2004) for recent reviews) were dis- 
covered by their recurrent emission of soft (<100 keV) 



Send offprint requests to: D. Gotz, email: diego@mi.iasf.cnr.it 
* Based on observations with INTEGRAL, an ESA project 
with instruments and science data centre funded by ESA 
member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain), Czech Republic and 
Poland, and with the participation of Russia and the USA. 



gamma-ray bursts. They undergo sporadic, unpredictable 
periods of activity, which last days to months, often fol- 
lowed by long periods (up to years or decades) during 
which no bursts arc detected. These recurrent bursts 
have typical durations of ^0.1 s and luminosities in 
the range 10 39 -10 42 erg s _1 . Occasionally, SGRs also 
emit giant flares that last up to a few hundred seconds 
and exhibit remarkable pulsations that reveal their spin 

Hurley et al. (1999)1 



periods (e.g. Mazets et al. (1979) 



Hurley et al. (2005) I 
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The spin period of SGRs can also be measured in their 
persistent (quiescent) X-ray emission. Typical luminosi- 
ties (0.5-10 keV) of these heavily absorbed sources are 
of the order of a few 10 35 erg s _1 . Persistent emission 
from SGR 1806-20 has been detected recently also in the 
soft 7-ray range up to ~150 keV ( |Mereghetti et al. 2005a| 
IMolkov et al. 2005|l . 

The bursting activity and the persistent emis- 
sion are generally explained in the framework of the 
"magnetar" model (see e.g.|Duncan fc Thompson (1992) 



Paczynski (1992)1 |Thompson fc Duncan (1995)| ), which 



involves slowly rotating (P ~ 5-8 s) highly magnetised 
(B ~10 15 G) isolated neutron stars. The magnetar model 
is based on the fact that the rotational energy loss in- 
ferred from the spindown (in the 10 -10 -10 -13 s s _1 range) 
is not sufficient to power the persistent X-ray luminosity. 
Hence it is the decay of the magnetic field itself that pro- 
vides the necessary energy. In this framework the mag- 
netic dissipation causes the neutron star crust to frac- 
ture. These fractures generate sudden shifts in the mag- 
netospheric footprints, which trigger the generation of 
Alfven waves, which in turn accelerate electrons above 
the pair-production threshold, resulting quickly in an opti- 
cally thick pair-photon plasma. The cooling of this plasma 
generates the typical short bursts of soft 7-ray radiation. 
An alternative explanation for short bursts, proposed by 



Lyutikov (2003) does not involve crustal fractures but 



simply local magnetic reconnection. The longer, much 
more energetic and rare flares are powered by a sudden 
reconfiguration of the magnetic field through magnetic re- 
connection (as in solar flares), and involve the entire neu- 
tron star magnetosphere. 

After a period of quiescence, SGR 1806-20 became 
active in the summer of 2003 Purley et al. 2003| ). 
Its activity then increased in 2004 (see e.g. 
IMereghetti et al. 2004a| IGolenetskii et al. 20"04|l . A 
strong outburst ~ ^Mereghetti et al. 2004b| ) during which 
about one hundred short bursts were emitted in a few min- 
utes occurred on October 5 2004 (see Section I3~2*|) . Finally 
a giant flare, whose energy (~10 46 erg) was two orders of 
magnitude larger than those of the previously recorded 
flares from SGR 0526-66 and SGR 1 900+14, was emitted 
on December 27 th 2004 (see e.g. |Palmer et al. (2005) 



Mereghetti et al. (2005b)[|Hurley et al. (2005)| )" 

Here we report on all the short bursts detected by 
INTEGRAL in 2003 and 2004. 

2. Observations and data analysis 

The observations of this source obtained with 
the INTEGRAL satellite in September and 
October 2003, during a period of moderate burst- 
ing activity, have been reported earlier (see 
Gotz et al. (2004a)| |Gotz et al. (2004b)|). We have 



re-analysed those data and, in addition, we report here 
on a much larger burst dataset. Excluding the bursts of 
the October 5 th event (see below), we have detected and 
analysed 224 bursts. The observations are summarised in 



Tab. 2] Our dataset represents the largest sample of short 
bursts imaged from SGR 1806-20 to date. As can be seen 
from Fig. the burst activity of SGR 1806-20 increased 
significantly with time, and the source was particularly 
active during Fall 2004. 
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Fig. 1. SGR 1806-20 burst rate vs. time. Only the bursts 
with fluence > 3x 10 -8 erg cm -2 are considered. The three 
time periods correspond roughly to Fall 2003, Spring 2004 
and Fall 2004. 



Here we concentrate on the data of the 
INTEGRAL imager IBIS <|Ubertini et al. 2003jl . We 
do not use the data of the X-ray monitor JEM-X 
l|Lund et al. 2 003). since only a small fraction of the 
bursts were inside its smaller field of view and, owing 
to their brightness, they were affected by saturation 
and dead time effects. All the bursts were in the field 
of view of the SPI instrument ( Ve drenne et al. 20 03). 
but at energies below ~300 keV SPI has a much smaller 
effective area compared to IBIS, which is therefore the 
most suitable instrument for detailed studies of short 
bursts. 

IBIS is a coded mask telescope with a large field of 
view (29°x29° zero sensitivity, 9°x9° full sensitivity). 
We used the low energy (15 keV-1 MeV) detector layer 
ISGRI IjLebrun et al. 2003|) . composed of 128x128 inde- 
pendent CdTe detectors (pixels) yielding a geometric area 
of ^1300 cm 2 on axis. ISGRI operates in a photon-by- 
photon mode: this means that for each photon the en- 
ergy, arrival time and interaction pixel are known. To 
identify the bursts we used the trigger information pro- 
vided by the INTEGRA L Burst Alert System (IB AS; 
Mereghetti et al. (2003) \ for the Target of Opportunity 
(ToO), public and Core Program data. For the Galactic 
Centre Field data, we searched for the bursts by comput- 
ing light curves with 10 ms time resolution and looking 
for significant excesses corresponding to the direction of 
SGR 1806-20 (for details see IMolkov et al. 2005p . All the 
bursts are unambiguously associated, within the typical 
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Excluding the huge outburst of October 5 (see text) 
Table 1. Observation summary of the bursts from SGR 1806-20 analysed here. (P: Public Data, CP: Core Program 
Observations, TOO: Target of Opportunity Observations, R: Open time Galactic Centre Data) The double horizontal 
line indicates the separation between the periods of different activity of the source (see text). 



~2' location uncertainties, with the position of the X-ray 
counterpart of SGR 1806-20 ( |Kaplan et al. 2002| ). 

For each burst we derived the 15-100 keV light curve 
with 10 ms time resolution. This energy band was chosen 
because no emission above 100 keV is detected in the typ- 
ical short bursts. In order to increase the signal-to-noise 
ratio, the light curves have been extracted by selecting 
only the pixels illuminated by the source over at least 
half their surface. The background was estimated by fit- 
ting a constant to 2 s time intervals before and after the 
burst. The light curves were then corrected for vignetting, 
caused by the fact that the bursts were detected at differ- 
ent off-axis angles, and for ISGRI dead time. The detector 
is made of 8 identical modules and the good events dead 
time, r, for each module is ~114 /is. For each burst, based 
on its position in the IBIS field of view, we determined the 
number of modules involved in the detection for at least 



half their surface, M, and applied the following dead time 
correction 

where N rec is the recorded flux (cts/s) in each time 
bin and N inc is the reconstructed incident flux. By evalu- 
ating this relation, one sees that only for the bins with a 
ratio between the recorded flux and the maximum record- 
able flux (M/r) larger than 0.3 a significant correction is 
required. For the rest of the bins the dead time effect is 
less than a few percent. It turns out that only 5.3% of the 
bursts have peak fluxes above this limit and their light 
curves have been corrected accordingly. 

These light curves have been used to determine the 
Tg durations (i.e. the time during which 90% of the burst 
fluence is accumulated), the peak fluxes (counts/bin) and 



4 



D. Gotz et al.: SGR 1806-20 bursts observed with INTEGRAL 



the fluences (counts) of each burst. To convert the peak 
fluxes and fluences to physical units we have used the con- 
version factor 1 count s _1 = 1.5xl0~ 10 erg cm -2 s _1 , 
derived by the spectral analysis of the brightest bursts 
(see Gotz et al. (2004a) I. These bursts are well fit with a 
Thermal Bremsstrahlung (TB) model with temperatures 
(kT) between 30 and 40 keV. The above conversion factor 
is valid for the average temperature of 32 keV. A varia- 
tion of the temperature between the two extremes implies 
a variation of only 7% of the conversion factor. So we have 
applied it also to fainter bursts, under the assumption that 
the spectral shape is the same for all the bursts. We point 
out, however, that extracting the average spectrum of the 
2003 bursts and comparing it to the average spectrum of a 
representative subsample of the bursts of Fall 2004, during 
which the source was more active, we find that the latter 
is marginally harder, but still within the above mentioned 
limits. Both average spectra are well represented by the 
TB model. 
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Fig. 2. IBIS/ISGRI light curves in the soft (20-40 keV, 
upper panel) and hard (40-100 keV, middle panel) en- 
ergy range and hardness ratio (lower panel) for a short 
burst from SGR 1806-20. Time=Q corresponds to 2004- 
08-18T21:52:49.39 UT 



3. Results 

3.1. Spectral Evolution 

Evidence for spectral evolution within the bursts of 
SGR 1806-20 was reported for the first time by 
Gotz et al. (2004a) who found that some bursts show 
significant spectral evolution, while others, particularly 
those with a "flat topped" profile, do not. In addition, 
a hardness-intensity anti-correlation within the bursts of 
SGR 1806-20 was reported by the same authors. 

We investigated the spectral evolution of the bright- 
est bursts by computing the background subtracted light 
curves in two different energy bands (20-40 (S) and 40-100 
(H) keV) and evaluating the hardness ratio, HR=(H- 
S)/(H+S), in time bins with the same number of total 
counts. The energy bands used here are slightly different 
from those used in Gotz et al. (2004a) because we realized 
that the individual pixels behave differently at energies be- 
tween 15 and 20 keV. In our analysis the spectral evolution 
is confirmed: excluding "flat topped bursts" , peaks tend to 
be softer than tails. Two good examples are shown in Figs. 
[21 and |3 However, a number of bursts show a less certain 
spectral evolution. By considering all the individual bins 
of all the bursts we derived the hardness (HR) intensity 
(/) distribution shown in Fig.0] To investigate if the data 
are correlated, we computed the Spearman rank-order cor- 
relation coefficient of the 217 data points, R s , which is - 
0.49. This corresponds to a chance probability of 4xl0 -15 
(7.4 a) that our distribution is due to uncorrelated data. 
According to an F-test the data are significantly (8a) bet- 
ter described by a linear fit (HR = 0.47 - 0.22 x log(J)) 
than by a constant. The parameters of the fit are very 
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similar to those found by Gotz et al. (2004a) So we can 
conclude that an anti-correlation between hardness and 
intensity is confirmed by our data. The inclusion or ex- 
clusion of flat-topped events from our analysis does not 
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Fig. 3. IBIS/ISGRI light curves in the soft (20-40 keV, 
upper panel) and hard (40-100 keV, middle panel) en- 
ergy range and hardness ratio (lower panel) for a short 
burst from SGR 1806-20. Time=Q corresponds to 2004- 
08-21T03:38:26.31 UT. 



affect our results, since they represent a small fraction of 
the bursts. 



3.2. The Large Outburst of October 5 2004 

On October 5 2004 IBAS triggered at 13:56:49 UT on the 
first of a series of bursts emitted from SGR 1806-20 . The 
source remained active, emitting several tens of bursts, 
until 14:08:03 UT. 

Fig. El (lower panel) displays the initial and brighter 
part of the ~ 11 minute long outburst. The time scale 
starts at 13:55:19 UT. We will refer to this time scale 
throughout our analysis. Two bright clusters of bursts are 
visible at t~100 s and t~280 s. They are so bright that the 
satellite telemetry was partially saturated (0 counts in the 
plot) for about 10 and 20 seconds respectively and only 
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Fig. 5. Light curves of the initial part of the October 5, 2004 outburst of SGR 1806-20. Upper panel: light curve at 
energy greater than ^80 keV obtained with the SPI Anti-Coincidence System in bins of 0.5 s. Bottom panel: light 
curve in the 15-200 keV energy range obtained with the IBIS/ISGRI instrument (bin size 0.1 s). The gaps in the 
IBIS/ISGRI light curve are due to saturation of the satellite telemetry. 
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Fig. 4. Hardness ratio ((H-S)/(H+S)) versus total count 
rate (20-100 keV, corrected for vignetting). The points 
are derived from the time resolved hardness ratios of the 
bursts with the best statistics. The line indicates the best 
fit with a linear function given in the text. 



part of the registered events could be sent to the ground. 
The two bright clusters are shown in more detail in Figs. 
[7|and|S] Despite the presence of the telemetry gaps in the 
data, one can see that they consist of many short bursts, 
with significant variability down to ~10 ms time scales. 
At least 75 short bursts, typical in terms of duration and 
shape, have been detected in IBIS/ISGRI data. 

We have extracted a spectrum over the entire duration 
of the outburst (676 s integration time, see Fig.EJ). A fit 
in the 20-300 keV band with a thermal bremsstrahlung 
model using XSPEC vll.2 IjArnaud et al. 1996|> gives a 
good result (\ 2 — 39.64 for 40 dof, allowing for 5% system- 
atics) and yields a temperature of 58±2 keV, higher than 
that usually measured for short bursts. A slightly worse 
but still acceptable fit (\ 2 = 45.99 for 40 dof) can be ob- 
tained using the sum of two black bodies, as proposed by 



Feroci et al. (2004) for the cumulative spectrum of short 



bursts from SGR 1900+14 and by |Qlive et al. (2004)| for 
an intermediate duration burst from the same source. The 
temperatures of the two black bodies are kTi=5.8±0.4 
keV and kT 2 = 18±0. 9 keV. Both valu es are higher 
than those derived by |Feroci et al. (2004)| (kTi=3.3±0 .1 
keV, kT 2 =9.5±0.4 keV) and by |Qlive et al. (200 
(kTi=4.3±0.1 keV, kT 2 =9.9±0.3 keV). On the other 
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Fig. 6. Spectrum of the October 5 outburst of SGR 
1806-20 measured with the IBIS/ISGRI instrument. 
Upper panel: data points and best fit with a thermal 
brcmsstrahlung model with kT=58 keV. Lower panel: 
residuals from the best fit in units of sigma. 



a mixture of different states, which may also explain the 
failure of the two black body model. 

We have extracted 5 more spectra for 2 bright bursts 
at t = 205.45 and t = 237 s and for 3 short (-0.1 s) 
bursts at t = 406.4, 560.25 and 761.8 s respectively. All 
the spectra are better represented by TB models than 
by power laws. The derived temperatures are reported 
in Tab. [3 Although no specific trend can be identified, 
there are significant variations between the bursts' tem- 
peratures. The spectra of the bursts at t = 205.45, 237 
and 560.25 s have very high temperatures, above 40 keV, 
which is the maximum temperature usually measured for 
short bursts. These events are similar in spectral hardness 
to two bursts detected from SGR 1900+14 following its 
giant flare H Woods et al. 1999|) . 
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0.2 


34.5±11 
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0.45 


65±7 


2.6xl0" 7 


23/20 


761.8 


0.2 


44±7 


2.8xl0 -7 
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The fluxes of the two bright clusters are not corrected for 
telemetry saturation. 

Table 2. Temperatures, fluxes and statistical goodness 
of the fits of the different bursts obtained with the TB 
model. 



hand by fixing the temperatures to the values derived by 
these authors we get unacceptable fits. This is possibly 
due to the fact that their spectra extend to lower energy. 
In the ISGRI energy range this model does not improve 
the reduced x 2 of the fit and hence the results of the two 
black body fit have to be taken with care. 

We computed the hardness ratios over the time pe- 
riods of the two clusters of bright bursts (t = 88-108 s 
and t = 260-293 s respectively). They indicate a clear 
hardening for the second cluster with respect to the first 
one. Even though the ISGRI detector partially saturates 
over these time intervals, the spectra still have enough 
counts to perform spectral analysis (see Figs. and 0. 
We have hence extracted two spectra over these two time 
periods. Applying TB models to the spectra, we see that 
the temperature increases significantly with time. In fact 
for the first cluster of bursts we derive a temperature (kT) 
of 39±2 keV, while for the second cluster we get 55.5±1.5 
keV. We point out however (see Tab.[5J) that the TB model 
does not fit the spectrum of the first cluster well. We di- 
vided this period into smaller time intervals and found 
that it is only the initial part (t = 89-92 s), that can- 
not be represented by a TB model. Its spectrum is in fact 
much better fitted (x 2 /dof = 36.6/31) by a power law with 
photon index r=2.80±0.05. The overall spectrum is hence 



The fluence of the entire outburst as measured by 
ISGRI is 1.5xl0 -5 erg cm -2 . This value is however heav- 
ily affected by the saturation of the brightest bursts and 
represents only a lower limit to the real fluence. In order to 
derive the overall fluence, we used the data provided by the 
Anti-Coincidence System (ACS) of the INTEGRAL spec- 
trometer (SPI). As can be seen in Fig. [5] (upper panel), 
only the brightest bursts are visible in these data and 
hence they represent complementary information to the 
ISGRI data. 

ACS data consist of the total count rate (50 ms time 
resolution) above 80 keV measured by the 91 bismuth 
germanate (BGO) scintillator crystals that surround the 
INTEGRAL spectrometer. The crystals are used as the 
anti-coincidence system of the spectrometer but are also 
capable of detecting high-energy transient events such as 
bright GRB and SG R bursts. We used the Mont e Carlo 
package MGGPOD ( |Weidenspointner et al. 2005| ) and a 
detailed mass modelling of SPI and the whole satellite 
dWeidenspointner et al. 2003| and references therein) to 
derive the effective area of the ACS for the direction of 
SGR 1806-20 . We computed the ACS light curve with a 
binsize of 0.5 s and estimated the background by fitting 
a constant value to all the data of the same pointing ex- 
cluding the bursts. We used the background subtracted 
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Fig. 7. IBIS/ISGRI light curves (15-200 keV) of the first cluster of bright bursts. Upper Panel: overall light curve with 
0.05 s resolution. The lower panels show some parts of the light curve in detail with a time resolution of 5 ms. The 
time scale is in seconds and is the same as for Fig. [5] 



light curve to compute the fluence of each burst cluster in 
counts. The ACS data do not provide any spectral infor- 
mation, so we computed the conversion factor to physical 
units based on the spectral shapes derived from ISGRI 
data and on the effective area computed through our simu- 
lations. The resulting fluences above 80 keV are 1.2xl0 -5 
and 9.4xl0 -6 erg cm~ 2 for the first and second clusters 
respectively. Converting these fluences to the 15-100 keV 
band one obtains 7.4xl0~ 5 and 3.2xl0~ 5 erg cm -2 re- 
spectively. By adding these results to the ones obtained 
for the ISGRI total spectrum, one can derive the total en- 
ergy output during the whole event, which is 1.2xl0~ 4 erg 
cm~ 2 . This corresponds to 3.25 xlO 42 erg for an assumed 
distance of 15 kpc ( |Corbel fc Eikenberry 2003| ). 



3.3. Number-Intensity Distribution 

The fluences of the 224 bursts derived above (excluding 

-th 



ent off-axis angles. The faintest bursts are missed when 
the source is observed at large off- axis angles. In order 
to correct for this effect we have computed the effective 
exposure of the source, taking into account the variation 
of sensitivity at various off-axis angles. This yields the 
exposure-corrected cumulative distribution shown by the 
dashed line in Fig. [5| 

Since the numbers at each flux level are not sta- 
tistically independent, one cannot use a simple % 2 
minimisation approach to fit the cumulative number- 
intensity distribution. So we have used the unbinned de- 
tections and applied the Maximum Likelihood method 
l|Crawford et al. 19 70) , assuming a single power-law dis- 
tribution for the number-flux relation (N(> S) oc S~ a ). 
We have used only the part of the distribution where com- 
pleteness was achieved (i.e. S >3xl0 -8 erg cm -2 ). In this 
case the expression to be maximised is 



the ones detected during the October 5 event) have been C = T\na — a In Si — T ln(l — 6~") 



used to compute the number-intensity distribution (Log 
N-Log S) of the bursts. 

The experimental distribution deviates significantly 
from a single power-law (Fig. |5J. This is first of all due 
to the fact that the source has been observed at differ- 



(2) 



where Si are the unbinned fluxes, b is the ratio between the 
maximum and minimum values of the fluxes, and T is the 
total number of bursts. This method yields a=0.91±0.09. 
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Fig. 8. IBIS/ISGRI light curves (15-200 keV) of the first cluster of bright bursts. Upper Panel: overall light curve with 
0.05 s resolution. The lower panels show some parts of the light curve in detail with a time resolution of 5 ms. The 
time scale is in seconds and is the same as for Fig. [5] 
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Fig. 9. Number-intensity distribution of all the bursts de- 
tected by INTEGRAL in 2003 and 2004. The continuous 
line represents the experimental data, while the dashed 
line represents the data corrected for the exposure. 

If a single power-law model is an adequate representation 
of data, the distribution of the quantities 

1 - S7 a 



l 



(3) 



should be uniform over the range (0,1). In our case, a 
Kolmogorov-Smirnov (K-S) test shows that a power law 
is an appropriate model, yielding a probability of 98.8% 
that the data are well described by our model. 

We then divided the bursts in two samples comprising 
51 and 173 bursts respectively. The division is based on 
the periods of different activity. The first period ends at 
orbit 181 and the second one starts at orbit 225 (see Tab. 
^) . The two slopes derived with the Maximum Likelihood 
method are a=0.9±0.2 for the low level activity period 
and a=0.88±0.11 for the high level one. The two slopes 
are statistically consistent with each other and a K-S test 
shows that the probability that the two distributions are 
drawn from the same parent distribution is 93%. Thus 
we conclude that the the relative fraction of bright and 
faint bursts is not influenced by the level of activity of the 
source. 

4. Discussion 

The bursts detected by INTEGRAL and presented here 
are common bursts as far as many aspects are concerned, 
such as durations (~0.12 s on average), peak fluxes, flu- 
ences and spectra but are the faintest detected at these 
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energies. Thanks to imaging we are confident that all of 
them were emitted by SGR 1806-20 . 

The good quality of our data has allowed us to study 
these bursts in detail. In particular, with respect to previ- 
ous experiments, we have a better combination of sen- 
sitivity, timing and spectral capabilities also for faint 
events. Hence we have been able to confirm the early find- 



ings (see Gotz et al. (2004a) I on the spectral evolution of 
SGR 1806-20 short bursts. The fact that burst peaks tend 
to be softer than their tails leads to an overall hardness- 
intensity anti-correlation. This characteristic still has no 
clear explanation within the magnetar model. 

The relatively large number of bursts has allowed us 
to constrain the shape and slope of the fluence distribu- 
tion, which is well described by a single power-law on 
over 2.5 orders of magnitudes. The INTEGRAL LogN- 
LogS has a slope of 0.91±0.01, similar to that derived 
at larger fluences (10~ 7 -10~ 5 erg cm -2 ) with KONUS 
data ( |Aptekar et al. 2001| ). This slope is also compati- 
ble, within 2 a, with the 0.71±0.11 obtained by combin- 
ing BATSE and ICE data flGogiis, et al. 2000| ) over the 
5xl0~ 8 -10 -5 erg cm~ 2 fluence range. On the other hand 
our LogN-LogS is significantly steeper than that obtained 
with RXTE, which has a slope of 0.43±0.06 over th e 
10~ 10 -10~ 7 erg cm~ 2 fluence range, (Gog iis, et al. 2 000). 
The difference between our slope and the RXTE one, 
which extends to lower fluences, is significant: this may 
imply that there is a break in the distribution, but we 
point out that in the overlapping part of the two datasets 
(6xl0^ 9 < 5" < 10~ 7 erg cm" 2 ) the RXTE slope is clearly 
statistically rejected by our data. We can therefore con- 
clude that, considering IBIS and BATSE//Ci?data, a sin- 
gle power-law gives a good representation of the cumula- 
tive energy distribution of SGR 1806-20 bursts over 3.5 
orders of magnitude. 

Persistent emission in the 20-150 keV range has been 
detected from SGR 1806-20 by integrating over long 
time intervals part of the IBIS/ISGRI data reported here 
HMereghetti et al. 2005a] IMolkov et al. 2005jl . Although 
the persistent emission has a spectrum harder than that 
of typical bursts, it is worth to check whether a significant 
part of such persistent emission, could be due to the cu- 
mulative contribution of numerous bursts too faint to be 
detected individually. Integrating our LogN-LogS distri- 
bution down to S~10 -10 erg cm -2 we find that the burst 
contribution is at most a few percent (~10~ 12 erg cm~ 2 
s _1 ) of the total. 

The hard and very energetic outburst of October 
5 th resembled a similar event seen in SGR 1900+14 
by KONUS on May 30 1998 < |Aptekar et al. 2001| . In 
that case the outburst lasted 5 minutes and a total of 



26 bursts were detected. The fluence was 3.5x10 4 



erg 



cm" , similar to what we measure here, and the temper- 
atures of the single bursts showed no particular correla- 
tion with time, but were on average at the high end of 
their distribution. In both sources, these energetic out- 
bursts formed by a rapid sequence of relatively hard short 
bursts, preceded of a few months the occurrence of gi- 



ant flares. The October 5 th event fits in the trend of in- 
creasing source activity shown by SGR 1806-20 in the 
last two years and also manifested in the rise in lumi- 
nosity and spectral hardness of the persistent emission at 
high (20-150 keV, |Mereghetti et al. 2005aD and low (2-10 
keV, |Mereghetti etaL^005c| ) energies. On the other hand, 
this peculiar event did not mark a peak or a turnover in 
the SGR activity. In fact the two XMM observations of 
SGR 1806-20 performed just before (September 6 2004) 
and the day after this large outburst (as a ToO in re- 
sponse to it) yielded similar spectral parameters, fluxes 
and pulse profiles, and bursts were seen in both observa- 
tions flMereghetti et al. 2005c| ). Also the INTEGRAL data 
indicate that the bursting activity remained high after this 
event. This also happened in 1998 to SGR 1900+14, which 
after the May 30 event remained in a very active state 
leading to the August 27 giant flare ( |Aptekar et al. 2001| ) . 
These results can be explained in the framework of a 



recent evolution of the magnetar model: Lyutikov (2003) 
explains SGR bursts as generated by loss of magnetic equi- 
librium in the magnetosphere, in close analogy to solar 
flares: new current-carrying magnetic flux tubes rise con- 
tinuously into the magnetosphere, driven by the deforma- 
tions of the neutron star crust. This in turn generates an 
increasingly complicated magnetic field structure, which 
at some point becomes unstable to resistive reconnection. 
During these reconnection events, some of the magnetic 
energy carried by the currents associated with the mag- 
netic flux tubes is dissipated. The large event described 
here can be explained by the simultaneous presence of 
different active regions (where the flux emergence is es- 
pecially active) in the magnetosphere of the neutron star. 
In fact, a long outburst with multiple components is ex- 
plained as the result of numerous avalanche-type recon- 
nection events, as reconnection at one point may trigger 
reconnection at other points. This explains the fact that 
the outburst seems to be composed by the sum of sev- 
eral short bursts. This kind of event might indicate a par- 
ticularly complicated phase of the magnetic field struc- 
ture which eventually led to a global restructuring of the 
whole magnetosphere with the emission of the giant flare 
on December 27. This mode also suggests that short events 
are due to reconnection, while longer events have in addi- 
tion a large contribution from the surface, heated by the 
precipitating particles, and are harder. This may explain 
the generally harder spectra observed. 

Thus events like these release a small (compared to 
giant flares) fraction of the energy stored in the twisted 
magnetic field of the neutron star, not allowing the mag- 
netic field to decay significantly. They are rather related to 
phases of high activity due to large crustal deformations 
(indicating that a large quantity of energy is still stored 
in magnetic form) and can be looked at as precursors of a 
major reconfiguration of the magnetic field. 

Finally we would like to point out that after the gi- 
ant flare of December 27 2004, SGR 1806-20 has re- 
mained active, with a burst rate per day of ~ 1.4 be- 
tween February 16 and April 28 2005. This prelimi- 
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nary value indicates that the level of activity of the 
source is intermediate between 2003 and 2004. The emis- 
sion of the giant flare, triggered by the reconfiguration 
of the magnetic field, has lowered the crustal stresses 
due to the magnetic dissipation ( |Thompson et al. 2002| ). 
A similar behaviour is seen in the X-ray band, 
where the persistent flux is intermediate between the 
2003 and the 2 004 observations <|Mereghetti et al. 2005c| 
IRea et al. 20051 |Tiengo et al. 2005| ).~ 

5. Conclusions 

We have presented the results of two years of 
INTEGRAL monitoring of SGR 1806-20 . During this 
time period the source went from a 7-ray quiescent state 
into a very active one that culminated in the giant flare 
of December 27, 2004. 

Our good quality data for low fluence bursts allowed 
us to establish that the number-density relation of the 
bursts is well represented by a single power law with 
index a=0.91±0.09. Despite the increase in the rate of 
emitted bursts, the burst properties do not change signif- 
icantly with time, neither does the slope of the number- 
density distribution. The spectral evolution of the bursts 
discovered in the 2003 data is confirmed also by this im- 
proved analysis of a larger sample of 224 bursts. The fact 
that bursts peaks tend to be softer than tails results in a 
hardness-intensity anti-correlation within the bursts. 

On October 5, 2004, close to the peak of its activity 
SGR 1806-20 emitted within 10 minutes more than a hun- 
dred short bursts that had a spectrum slightly harder than 
usual and involved a total energy release of 3xl0 42 erg. 
These high temperature bursts can be explained in the 
framework of the magnetar model as an avalanche- type 
reconnection event in the neutron star magnetosphere 
caused by a particularly complicated structure of the mag- 
netic field. Events of this kind might be precursors of the 
major reconfiguration of the whole magnetic field of the 
neutron star causing large flares in SGRs. 
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